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Abstract 

The measurement of D meson production provides key tests for parton energy-loss models, which 
predict that charm quarks should experience less in-medium energy loss than light quarks and 
gluons. The ALICE experiment has measured the production of prompt D", and D*+ mesons 
in pp and Pb-Pb collisions at the LHC at = 7 and 2.76 TeV and at -^Snn = 2.76 TeV, respec- 
tively, via the exclusive reconstruction of their hadronic decay. The /7T-differential production 
yields in the range 2 < pj < \6 GeV/c at central rapidity, \y\ < 0.5, were used to calculate the 
nuclear modification factor. A suppression of a factor 3 to 4 for transverse momenta larger than 
5 GeV/c in the 20% most central collisions was observed. Preliminary results in an extended 
px-range, using the data sample collected during the 201 1 Pb-Pb run, together with the first 
measurement of nuclear modification factor will be shown. 



1. Introduction 

Under conditions of high-energy density and temperature, produced in ultra-relativistic heavy- 
ion collisions, lattice QCD calculations predict the transition from ordinary nuclear matter to a 
deconfined state of quarks and gluons called Quark Gluon Plasma (QGP) [1!- Heavy-flavour 
quarks, differently from light quarks and gluons, are only produced at the early stage of the col- 
lision in high-virtuality scattering processes. They traverse the medium and are expected to be 
sensitive to its density through the mechanism of in-medium partonic energy loss. Heavy flavours 
should lose less energy than light-quarks and gluons as a consequence of a mass-dependent re- 
striction in the phase space into which gluon radiation can occur f2\. Moreover, if in-medium 
hadronization is the dominant mechanism of charm-hadron formation then the strange charmed 
hadrons, like D^, should be largely enhanced Ol. The nuclear modification factor (7?aa) of D 
mesons, obtained by comparing their production in proton-proton and heavy ion collisions, al- 
lows to probe the aforementioned properties. The /?aa is defined as: R^j^ifj) — -^jt^ dcrD^/dp^ 

where A^^^ is the yield in A-A collisions, (Taa), in a given centrality class, is the average nu- 
clear overlap function calculated via Glauber model and cr^^ is the production cross section of 
D mesons in pp collisions. In absence of medium effects is expected to be I. Using the 
data collected in 2010 we published the first measurement of the suppression of high transverse 
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momentum D mesons in central Pb-Pb collisions at -\/snn - 2.76 TeV (4]. We observed a sup- 
pression of a factor 3 to 4 for transverse momenta larger than 5 GeV/c in the 20% most central 
collisions. In this contribution we present the nuclear modification factor of D", D^, D*^ and 
as function of transverse momentum in the centrality range 0-7.5% as extracted from the data 
collected in 201 1. A comparison with model predictions is made. 

2. D meson reconstruction with ALICE 

The production of D°, D+, D*+ HISl and EEl was measured in pp and Pb-Pb collisions 
at central rapidity {\y\ < 0.5) via the exclusive reconstruction of the decays D" — > K^tt^, — > 
K-7T+7T^, D*+ -> D";r+ -> K-jt+tt^ and D+ ^ (fm^ K-K+7r+. The analysis strategy for the 
extraction of the signal out of a large combinatorial background is based on the reconstruction of 
the D meson decay topology (D" in case of the D*^). In order to resolve the vertex of the decay 
the Inner Tracking System (ITS) fW\ provides the required resolution of a few tens of microns on 
the track position at the primary vertex. D meson candidates are defined from displaced tracks 
and selected by means of topological cuts. To further suppress the combinatorial background, 
particle identification (PID) on the D meson decay tracks is employed. PID is performed using 
the information on specific energy deposit in the Time Projection Chamber (TPC) and on the 
time of flight measured by the Time of Flight (TOF) detector It allows a reliable identification 
of kaons and pions in a wide momentum range up to 2 GeV/c. The signal yield is extracted 
by fitting the invariant mass distribution using a Gaussian function for the signal peak. The 
background is fitted using an exponential shape in the case of D*\ and while in the case 
of the D*^ a threshold function convoluted with and exponential is chosen. The correction for 
efficiency and acceptance is performed using Monte Carlo simulations based on Pythia Perugia-0 
tuning and HIJING event generator. In order to extract the prompt D fraction, the contribution 
of D mesons from B decays was evaluated relying on FONLL [9 1 calculations, which describe 
well B hadron production at Tevatron |10| and LHC lfm[T2l . A data driven method, based on 
the different shape of the impact parameter to the primary vertex for primary and secondary D 
mesons is under study. The analyses presented in this paper are based on the Pb-Pb data sample 
at centre-of-mass energy y^NN = 2.76 TeV collected in November 2011. The events were 
triggered with centrality-based selection using information from Silicon Pixel Detector {\ri\ < 2) 
and the VZERO |8| scintillator hodoscopes (2.8 < ?; < 5.1 and -3.7 < rj < -1.7). Only events 
with a vertex found within 10 cm from the centre of the detector along the beam line were used, 
for a total of 16 x 10^ collisions in the 0-7.5% centrality class. 

3. D meson pp reference 

The D-meson production cross section in pp collisions at = 2.76 TeV, needed for the 
Raa calculation, is based on a scaling of the measured D meson production cross section at 
yfs - 1 TeV using the ratio of FONLL predictions at V* = 2.76 and 7 TeV. The procedure is 
validated with data by comparing the scaled cross section with the one measured on a pp sample 
with limited statistics at = 2.76 TeV. A detailed explanation of the procedure can be found 
in m. At high pj, 24 < pj < 36 GeV/c for D+, D*+ and 16 < pj < 24 GeV/c for D° a 
measured reference is not available. Therefore the cross section was extrapolated on the basis of 
the FONLL/data ratio observed at lower transverse momentum. In 1 < < 2 GeV/c, where 
the total uncertainties of the 7 TeV scaling and of the 2.76 TeV measurement are comparable, 
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the reference was calculated as the weighted average of the 7 TeV scaled and of the 2.76 TeV 
measurement using the relative uncorrelated uncertainties as weights. 

4. D meson nuclear modification factors 

The nuclear modification factor of D°, D^, D*^ and mesons measured with the data taken 
in 2011 in the centrality class 0-7.5% is shown on the left panel of Fig. [T] The of the four 
D mesons agrees within uncertainties and shows a clear suppression of about a factor 5 for 
Pj ~ 10 GeV/c. The estimated total systematic uncertainty, accounting for the uncertainties on 
signal extraction procedure, PID selection strategy, track reconstruction efficiency, cut stability 
and hypothesis on B mesons Raa (accounting for all the potential nuclear and medium effects 
affecting B production) varies from ~70% to 12% depending on p-Y and D meson species. The 
/?AA0fD°, D+, andD*+ were averaged. The contributions of the three D mesons to the average 
were wighted by their statistical uncertainties. The systematic errors were calculated by propa- 
gating the uncertainties through the weighted average, where the contributions from the tracking 
efficiency, B feed-down correction and scaling of the pp reference were taken as fully corre- 
lated among the three D meson species. On the right panel of Fig. [T] we compare our result with 
charged pions and charged hadrons Raa as measured by ALICE with 2011 data sample in the 
centrality class 0-10%. The three measurements are compatible within uncertainties, however 
a hint of a lower suppression of the D mesons is seen at low pj. The data points show a hint 
for R^j^ < 1 at low pj while at high pj the data are consistent with both a flat and a rising 
shape. Fig.|2]shows the average 7?^^ compared with models. We found that radiative energy loss 
supplemented with in-medium D meson dissociation |14| and radiative plus collisional energy 
loss in the WHDG 1 15 1 implementations, together with BDMPS-ASW 1 16 1 and POWLANG fnl 
describe the data reasonably well. The heavy quark transport model with in-medium resonance 
scattering and coalescence |18 | tends to underestimate the suppression observed in the data. 
We verified, using pQCD calculations based on the MNR |19| code and nuclear- modified PDF 
from the EPS09 |20| parametrization, that nuclear shadowing has a relatively small effect for 
Pj ~ 10 GeV/c. This indicates that the suppression, of about factor 5 at pj ~ 10 GeV/c indicate 
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Figure 1: Left panel: D", D+, D*+ and D+ meson nuclear modification factor in the centrality class 0-7.5% in the pj 
range 1 < px < 36 GeV/c (4 < px < 12 GeV/c in the case of DjJ^). Right panel: Comparison of the average D", and 
D** Raa with charged pions and charged hadrons. 
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Figure 2: Average Saa of D", and D*^ compared with model calculations 



a large energy loss of charm quarks in the hot and dense QCD medium formed in heavy-ion 
collisions at the LHC. 

5. Summary 

We have presented the first measurement of the D", D^, D*^ nuclear modification factor in 
the centrality class 0-7.5% using the 2011 Pb-Pb data. The pj range of the measurement is 
extended with respect to the 2010 result f4l and now covers the transverse momentum range 
1 < Px < 36 GeV/c. The p-Pb run scheduled for the beginning of 2013 will allow for the 
study of the initial state effect at low px. Finally we presented the first R^a of the meson, in 
the pt range A < pi < 12 GeV/c, together with its comparison with the other D mesons. The 
result of the comparison, not conclusive due to the large statistical and systematic uncertainties, 
is intriguing in view of future LHC runs. 
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